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ABSTRACT 

We report coincident damped Lyman alpha (DLA) and sub-DLA absorption at 
z a bs = 2.66 and z a b s = 2.94 towards the z ~ 3, 13.8 arcsecond separation binary 
quasar SDSS 1116+4118 AB. At the redshifts of the absorbers, this angular separation 
corresponds to a proper transverse separation of ~ 110 kpc. A third absorber, a 
sub-DLA at z a b s = 2.47, is detected towards SDSS 1116+4118 B, but no corresponding- 
high column density absorber is present towards SDSS 1116+4118 A. We use high 
resolution galaxy simulations and a clustering analysis to interpret the coincident 
absorption and its implications for galaxy structure at z ~3. We conclude that the 
common absorption in the two lines of sight is unlikely to arise from a single galaxy, 
or a galaxy plus satellite system, and is more feasibly explained by a group of two 
or more galaxies with separations ~ 100 kpc. The impact of these findings on single 
line of sight observations is also discussed; we show that abundances of DLAs may be 
affected by up to a few tenths of a dex by line of sight DLA blending. From a Keck ESI 
spectrum of the two quasars, we measure metal column densities for all five absorbers 
and determine abundances for the three absorbers with log N(HI) > 20. For the two 
highest N(HI) absorbers, we determine high levels of metal enrichment, corresponding 
to 1/3 and 1/5 Z@. These metallicities are amongst the highest measured for DLAs 
at any redshift and are consistent with values measured in Lyman break galaxies at 
2 < z < 3. For the DLA at z a b s = 2.94 we also infer an approximately solar ratio 
of a-to-Fe peak elements from [S/Zn] = +0.05, and measure an upper limit for the 
molecular fraction in this particular line of sight of log /(H2) < —5.5. 
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1 INTRODUCTION 

For the last decade, a simple schematic view of the rela- 
tive sizes of QSO absorption systems has been built on the 
coupling of absorbers with galaxies of varying luminosities 
and impact parameters and the consideration of luminosity 
functions and number densities (e.g. Steidel 1993; Steidel 
1995; Lanzetta 1993). Although a galaxy's gas cross-section 
depends on its individual properties (such as mass and lu- 
minosity), the topology of a given galaxy is usually consid- 
ered to be hierarchical and based on column density (e.g. 
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Steidel 1993; Churchill, Kacprzak & Steidel 2005). In this 
picture, the damped Lyman alpha (DLA) systems represent 
a relatively small fraction of a galaxy's gas cross-section, 
associated with its inner ~ 10-20 kpc. The Mg II bearing 
gas, associated with Lyman limit systems, occupies a some- 
what larger halo; in the original Steidel (1995) picture this 
halo was roughly spherical, had a covering factor of approx- 
imately unity and a radius ~ 40 kpc for an L* galaxy. The 
largest absorption cross-section was associated with C IV 
absorbing gas and extended out to distances on the order of 
100 kpc. Although attractive in its simplicity, this picture 
has recently undergone significant re-evaluation. For exam- 
ple, although the idea of large C IV halos has been vin- 
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dicated by observations of absorbers near to Lyman break 
galaxies (LBGs, Adelberger et al. 2003, 2005b), it has been 
argued that some CIV absorbers may be associated with a 
more diffuse component, possibly the intergalactic medium 
(Pieri, Schaye & Aguirre 2006). Structure has also been in- 
ferred in the Mg II population. Ellison et al. (2004) argued 
that variations in Mg II equivalent widths (EWs) on kpc- 
scales seen in spatially resolved lensed QSO images suggest 
individual Mg II 'clouds' that are an order of magnitude 
smaller than the halo sizes found by previous galaxy sur- 
veys. Even the sizes of the Mg II halos are currently being 
re-assessed; it has been suggested by Churchill, Kacprzak & 
Steidel (2005) that the original survey strategies may have 
led to an under-estimate of the extent of the absorbing gas. 

In light of these recent re-evaluations, it may be surpris- 
ing how little progress we have made in determining the sizes 
of DLAs. These are probably the best studied of the quasar 
absorption line menagerie, partially because of the numer- 
ous chemical elements that can be used to infer their star 
formation histories (e.g. Dessauges-Zavadsky et al. 2004). 
Unlike the case for C IV and Mg II absorbers, there is scant 
data from which we can directly infer DLA sizes. The two 
main techniques that have previously been used for this es- 
timation, namely the association of individual galaxies with 
absorbers and the application of lensed QSOs, are not yet on 
a solid statistical footing. The former of these techniques re- 
quires a relatively large sample of DLAs with known galaxy 
counterparts. Although Chen & Lanzetta (2003) attempted 
this with a sample of 6 galaxies, the statistics are poor and 
only available for z < 1, whereas the vast majority of known 
DLAs are at z > 2 (e.g. Prochaska, Herbert-Fort & Wolfe 
2005). At higher redshifts, only M0ller, Fynbo & Fall (2004) 
and Weatherley et al. (2005) have made direct detections of 
DLAs and find impact parameters in the range 2-25 ft™ 
kpc for 3 DLAs. Measurements of DLA sizes from lensed 
QSOs are currently limited by the very small number (4) of 
DLAs that have been detected in lensed sightlines and by the 
small transverse scales that they probe. Two out of the four 
cases (Churchill et al. 2003; Kobayashi et al. 2002) probe 
very small scales (< 250 pc), leaving only two measurements 
on kpc scales (d ~ 5 kpc by Lopez et al. 2005 and d ~ 10 kpc 
by Smette et al. 1995). A few indirect limits of DLA sizes 
also exist, for example lower limits based on extended back- 
ground radio emission (e.g. Foltz et al. 1988; Briggs et al. 
1989) or on the unique case of transverse Lya fluorescence 
discovered by Adelberger et al. (2006). It is also possible to 
estimate gas cross-sections by combining the DLA number 
density with a Holmberg relation between the radius of a 
disk and the galaxy luminosity (e.g. Fynbo et al. 1999). The 
disadvantage of this approach is that it assumes a specific 
and fixed geometry, although like other methods, it gives 
typical disk sizes up to ~ 30 kpc. Apart from this handful of 
direct and indirect constraints, the scenario in which DLAs 
have cross-sections of a few tens of kpc has been largely 
untested. 

In this paper, we present observations of a close binary 
QSO, SDSS 1116+4118 AB, hereafter QSO A/B, with an 
angular separation between the two components (both at 
z ~ 3) of 13.8 arcseconds, see Table [T] and Figure [T] QSO 
B exhibits three high column density, intervening absorbers, 
two of which are also detected in QSO A, even though the 
proper separation at the redshift of the absorbers is more 



than 100 kpc. In order to interpret this result, we use 
high resolution smoothed particle hydrodynamic (SPH) sim- 
ulations of two galaxies at z = 2.3, 3.0, 3.6 to assess whether 
a structure as large as 100 h? kpc is likely to be a single 
galaxy, a satellite system or a group of galaxies and discuss 
the implications of these possibilities. 

We adopt a cosmology of = 0.7, Q.m = 0.3 and Ho = 
70 km/s/Mpc. In this cosmology, redshifts of 2.47, 2.66 and 
2.94 (the redshifts of the absorbers studied in this paper) 1 
arcsecond in the transverse direction corresponds to proper 
linear distances of 8.09, 7.96 and 7.75 Kf kpc respectively. 



2 OBSERVATIONS AND DATA ANALYSIS 

Although the Sloan Digital Sky Survey (SDSS) spectro- 
scopic quasar survey has provided the largest sample (~ 10 ) 
of quasars in existence, it selects against close quasar pairs 
due to the finite size of optical fibres in the multi-object 
spectrograph. This fibre collision limit implies that only one 
member of a pair with AO < 1' will make it into the quasar 
catalog. Hennawi et al. (2006a) selected a large sample of 
candidate close companions around the SDSS quasars us- 
ing photometric redshift techniques, and spectroscopically 
confirmed them to be quasar pairs with follow-up observa- 
tions at low spectral resolution. A small number of QSO 
pairs can also be discovered when there is plate overlap 
for a given field. In this case, fibres can be placed on very 
closely spaced objects in two different plates. The QSO pair 
SDSS 1116+4118 AB falls into this latter category; an SDSS 
postage stamp of the field is shown in Figure [1] and shows 
both QSOs as well as a foreground galaxy. For a sample 
of confirmed spectroscopically confirmed QSO pairs, we are 
currently undertaking higher resolution spectra in order to 
study the transverse absorption properties of intervening 
galaxies and the intergalactic medium. The data presented 
here were obtained as part of that systematic study and they 
represent the only pair for which we have currently detected 
a DLA in either line of sight (many of the spectra do not 
have large Lya forest coverage). 

Sanchez- Alvaro & Rodriguez-Calonge (2007) have re- 
cently proposed SDSS 1116+4118 AB to be a wide separa- 
tion lens, with the foreground galaxy (whose photometric 
redshift is given as z = 0.25) as the lensing mass. However, 
we consider this unlikely due to a) the improved redshift 
determinations of the two QSOs (Table [TJ which are sig- 
nificantly different: the CIV emission lines are offset by 36 
A; b) the considerable differences between QSO A and B in 
our ESI spectra, both in absorption and emission character- 
istics; c) the extremely low likelihood that a single galaxy 
would produce such a wide separation image (the implied 
mass-to-light ratio is >100 at the Einstein radius); d) lack 
of counter images on north side of the lens galaxy (the off- 
set lens should produce a quadruple image); e) high redshift 
(z > 2.5) C IV absorbers show large fractional equivalent 
width differences (>> 50%) between QSO A and B, incon- 
sistent with lensing-predicted line of sight separations of ^ 
1.5 h^o kpc (e.g. Lopez et al. 2000; Rauch et al. 2001a,b; 
Ellison et al. 2004) and f) the inconsistency of QSO colours, 
e.g. g - i=0.61±0.02 and 0.45+0.02 for A and B respec- 
tively. Only the presence of a massive cluster would cause 
such wide separation multiple images, such as the recently 



Damped Lyman alpha absorption in a binary QSO 3 



i i i i i i i i i + 



♦ 



Figure 1. SDSS image of the binary QSO SDSS1116+4118 AB. 
QSO B is at the centre of the image and QSO A is offset by 13.8 
arcseconds to the north-east. The galaxy located between the two 
QSO images does not have a spectroscopic redshift, but Sanchcz- 
Alvaro & Rodriguez-Calonge (2007) give a photometric redshift 
of z = 0.25. 
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discovered 14 arcsecond separation quadruply imaged SDSS 
1004+4112 (Inada et al. 2003) where the estimated cluster 
mass is M ~ 10 14 M© (Oguri et al. 2004). No such clus- 
ter is seen in the SDSS images of SDSS 1116+4118, or in 
deeper MMT images (A. Marble, private communication). 
We therefore assume that SDSS 1116+4118 AB is a pro- 
jected pair of QSOs. 



2.1 Data Acquisition and Reduction 

On March 3 2006 we obtained 3300 seconds of integration 
in two exposures on QSO A/B using the Echellette Spec- 
trograph and Imager (ESI, Sheinis et al. 2002) on the Keck 
telescope. Since the entrance slit of the spectrograph is 20 
arcseconds long, the position angle was chosen so that both 
QSOs A and B were covered (see Figure [I}. The observing 
conditions were relatively poor, with high humidity and see- 
ing typically 1.3 arcseconds. The data quality were further 
degraded by broad absorption patterns across many echelle 
orders due to condensation on the dewar window caused by 
observing in high humidity. Nonetheless, with a 1 arcsecond 
slit width and lxl binning the S/N ratios per pixel were 
~ 50 in QSO A and 30 in QSO B at 6000 A. 

The data were reduced using a customized version of 
ESIRedujQ which was adapted to deal with multiple ob- 
jects on the slit. Extracted spectra were calibrated to a vac- 
uum heliocentric wavelength scale and determined to have 
FWHM resolutions of ~ 60 km/s (R ~ 5000). The spec- 
tra from the two individual exposures were combined by 
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Figure 2. Lya fit (dashed profile) and error (dotted profiles) to 
the z a bs ~ 2.47 sub-DLA seen towards QSO B (bottom panel), 
but not towards QSO A (top panel). The transverse separation 
of the two sightlines at this redshift is 112 kpc. 

weighting according to S/N. The QSO continuum (includ- 
ing the broad absorption features induced by the residue on 
the CCD) was estimated with the Starlink software DipscQ 
by fitting a cubic spline polynomial through unabsorbed re- 
gions of flux. To test the effect of the CCD residue, we mea- 
sure the equivalent widths of unsaturated absorption lines 
that fall in this region of the spectrum. In theory, dividing 
the spectrum by the continuum function should not alter the 
equivalents widths, if the fitting is accurate. We find equiv- 
alent widths that agree to within less than a few percent, 
indicating that the broad absorption does not affect our line 
measurements. 

We identified three high N(HI) absorbers present in at 
least one of the two lines of sight at redshifts z a b s = 2.47 (B 
only), 2.66 (A and B) and 2.94 (A and B). 

2.2 Column Density Determination 

The HI column densities of the absorbers at z&hs = 2.47, 
2.66 and 2.94 were determined by fitting the continuum 
normalized spectra with damped profiles using the Dipso 
software. The fits to the data are shown in Figures [2] [3] 
and [4] For all 5 absorbers, we detect metal lines, including 
Si II, Fe II, S II, Zn II and Ni II, see Figures E] © and 
Since the resolution of ESI is significantly larger than the 
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Table 1. Target Properties 



QSO 


RA (J2000) 


Dec (J2000) 


i mag 


•2cm 


SDSS 1116+4118 A 


11 16 11.7 


41 18 21.5 


17.97 


2.982+0.007 


SDSS 1116+4118 B 


11 16 10.7 


41 18 14.4 


19.00 


3.007+0.007 



1.5 
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Figure 3. Lya fit (dashed profile) and error (dotted profiles) to 
the z a b s ~ 2.66 DLA/sub-DLA seen towards QSO A (top panel) 
and QSO B (bottom panel). The transverse separation of the two 
sightlines at this redshift is 110 h^ kpc. 



Figure 4. Lya fit (dashed profile) and error (dotted profiles) to 
the z a t, s ~ 2.94 DLA/sub-DLA seen towards QSO A (top panel) 
and QSO B (bottom panel). The transverse separation of the two 
sightlines at this redshift is 107 h^ kpc. 



Doppler widths of the metal lines, Voigt component fitting 
does not yield a physically meaningful decomposition of the 
line profile. Instead, column densities are determined using 
the apparent optical depth method (AODM, e.g. Savage & 
Sembach 1991). Accurate column densities of typically weak 
lines such as Zn II can be well determined in lower resolu- 
tion spectra. For example, the early work of Pettini et al. 
(1997) determined N(Zn II) from spectra with FWHM ~ 
1 A, similar to ESI's resolution, which have subsequently 
been confirmed with higher resolution echelle spectroscopy. 
The AODM also allows us to assess the impact of saturation 
in unresolved lines when multiple transitions from a given 
species are observed. This is often the case for Fe II lines, 
and when saturation is suspected, we adopt the maximum 
column density, usually derived from the weakest transition. 
The main disadvantage of the AODM is that the assessment 
of contamination and blending is less straightforward. 

Due to their ionization potentials, it is usually assumed 
that the singly ionized species represent the dominant state 
of elements in DLAs. For absorbers with log N(HI) < 20.0 



there may be a non-negligible ionization correction which 
undermines this assumption. When a range of metal lines 
from different ionization states are detected, it is possible to 
model their relative contributions. Often this can be done 
with Al II plus Al III or Fe II plus Fe III. However, as de- 
scribed below, for the two absorbers in our spectra with log 
N(HI) < 20.0 such an estimate of ionization correction is 
not possible due to the limited number of transitions that 
we detect. In our abundance determinations we therefore as- 
sume that N(X) = N(X II) when log N(HI) > 20.0, but for 
the lower N(HI) absorbers we quote only column densities 
and not abundances. 

Fits to both HI and metal line species are now described 
in the following sub-sections on an absorber-by-absorber ba- 
sis. The redshifts that we allocate to each absorber refer to 
zero velocity for the metal lines, not the value determined 
from the Lya. 
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Table 2. Ionic Column Densities for z a t, s = 2.4675 sub-DLA towards QSO B 



Ion 


A 


/-value 


log N 


log N adopt 


H I 


1215.6701 


0.41640 




19.85±0.15 


Al II 


1670.7874 


1.880 


12.73±0.09 


12.73±0.09 


Fe II 


1608.4511 


0.0580 


<13.6 


13.33±0.12 


Fe II 


2260.7805 


0.00244 


<14.2 




Fe II 


2249.8768 


0.001821 


<14.4 




Fe II 


2344.2140 


0.1140 


13.33±0.12 




Fe II 


2374.4612 


0.0313 


14.10±0.09 




Fe II 


2382.7650 


0.3200 


13.42±0.07 




Si II 


1808.0130 


0.002186 


<14.4 


<14.4 


Zn II 


2026.1360 


0.4890 


12.08±0.26 


<12.3 


Zn II 


2062.6640 


0.2560 


<12.3 




Cr II 


2056.2539 


0.1050 


<12.6 


<12.6 
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Figure 5. Metal lines for the sub-DLA detected towards QSO B 
at 2 abs ~ 2.4677. 

2.2.1 z abs = 2.4675 sub-DLA towards QSO B 



Figure 6. A selection of detected metals for the DLA/sub-DLA 
detected towards QSO A (left column) and QSO B (right column) 
at 2 abs ~ 2.66. The x-axis shows a velocity scale relative to the 
centre of the metal line absorption, whose redshift is given in the 
bottom panels. 



Despite the low N(HI) of this absorber and the relatively 
poor S/N in the blue part of QSO B's spectrum, clear damp- 
ing wings allow us to constrain the hydrogen column density 
to within 0.15 dex (not including continuum fit errors). No 
corresponding DLA or sub-DLA is seen towards QSO A, al- 
though there is a saturated Lya absorber at z a b s = 2.4688 
corresponding to a velocity offset of Av — 95 km/s from the 
best fit HI redshift of z abs = 2.4677. 

Although the low N(HI) and relatively poor S/N pre- 
clude the detection of many metal line species, we do detect 



several lines of Fe II (although many are blended) as well as 
Al II A 1670 and Zn II A 2026, see Table We adopt the 
Fe II column density from the 2344 A line since the other 
detected transitions of this species (Fe II AA2374, 2382) are 
blended (see Figure [5). 

One potentially serious problem with lower resolution 
data is blending with either species at the same wavelength 
(e.g. Zn II A 2062 with Cr II A 2062) or from lines at other 
redshifts. At high resolution, this is usually easily identified 
during the Voigt profile fitting process, but is not always 
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Figure 7. A selection of detected metals for the DLA/sub-DLA 
detected towards QSO A (left column) and QSO B (right column) 
at z a bs ~ 2.94. The x-axis shows a velocity scale relative to the 
centre of the metal line absorption, whose redshift is given in the 
bottom panels. 



obvious in lower dispersion data. This issue is relevent in 
our determination of a Zn II column density, since the Zn II 
A 2026 line (the stronger of the two UV lines) is only sepa- 
rated from Mg I A 2026 by 50 km/s. The relative strengths 
of Mg I and Zn II will depend on a variety of factors, in- 
cluding intrinsic abundance ratios, differential dust deple- 
tion and ionization structure of the absorber. For DLAs and 
sub-DLAs with high column densities of metals, it seems 
that Mg I can contribute significantly to the total equiva- 
lent width at A ~ 2026 (e.g. York et al. 2006; Peroux et al. 
2006; Herbert-Fort et al. 2006). Although this is unlikely to 
be a serious problem for this sub-DLA, whose metal column 
densities are not large, we conservatively quote an upper 
limit for N(Zn II) based on the non-detection of Zn II A 
2062. Without a correction for Mg I and without an ion- 
ization correction (see below) the Zn abundance determined 
from Zn II A 2026 is high, about 0.4 Zq. It would therefore 
be very interesting to improve the S/N for this absorber to 
confirm (or otherwise) this relatively high metallicity. With 
the limited metal column densities we are able to reliably 
measure, it is not possible to determine ionization correc- 
tions for this absorber, so we do not attempt to convert our 
AODM measurements into abundances. 

It is interesting to note that although the extent of C IV 
galaxy halos are usually considered to have large sizes (up 
to a few hundred kpc, Steidel et al. 1995; Adelberger et al. 
2005b), in this case, there is a C IV absorber in QSO A at 



Zabs = 2.4532 compared with a redshift for C IV in QSO 
B (where the sub-DLA is detected) of z abs = 2.4683. This 
corresponds to a velocity difference of ~ 1300 km/s. It is not 
clear whether the C IV absorption in QSO A is associated 
with the sub-DLA detected in QSO B. The interpretation 
of C IV halo sizes clearly depends on the velocity tolerance 
permitted for line matches. We will analyse the C IV ab- 
sorbers in this pair, and a larger sample, in a forthcoming 
paper. 



2.2.2 z abs = 2.6617 DL A towards QSO A and z aba = 
2.6625 sub-DLA towards QSO B 



For the DLA towards QSO A we detect metal lines from 
the following species: Fe II, Si II, Zn II, Al III and Ni II and 
determine an upper limit for Cr II, see Table [3] and Figure 
[U] We also detect Al II A 1670, but this transition is blended 
with C IV at Zabs ~ 2.94 so we do not quote a column density 
for it. For the Fe II column density, we adopt an average of 
the AA 1608 and 2374 transitions which are the least likely 
to suffer from saturation. 

As mentioned in the previous subsection, Mg I may be a 
significant contaminant at 2026 A. The most common way to 
calculate the contribution from Mg I is to measure N(Mg I) 
from the Mg I A 2852 line and predict the contribution at 
2026 A according to the relative oscillator strengths. How- 
ever, at a redshift of z a b s ~ 2.66 the Mg I A 2852 line is 
shifted into the infra-red. We therefore rely on Zn II A 2062 
to provide a Zn II column density. Although the Zn II A 2062 
A line can itself suffer from blending with Cr II A 2062, in 
this case the upper limit on the stronger Cr II A 2056 rules 
out this possibility. Comparing the column densities derived 
from the two Zn II lines in Table [3] shows that the contribu- 
tion from Mg I at A 2026 would have led to an over-estimate 
of N(Zn II) by ~ 0.2 dex. 

For the sub-DLA towards QSO B, we only determine a 
column density for Fe II (see Table [4} since Al III AA 1854, 
1862 have very broad profiles and are slightly offset from 
the central velocity, so may be blends or mis-identifications. 
The usually strong Al II A 1670 is also in a blended part of 
the spectrum. 



2.2.3 z abs = 2.9^22 sub-DLA towards QSO A and z a bs = 
2.9467 sub-DLA towards QSO B 



The best fit N(HI) to the z abs ~ 2.94 absorber in QSO 
A yields a column density just below the DLA threshold, 
although within the error bars it may still be a 'classical' 
DLA. The absorber towards QSO B appears to be a blend, 
but the red damping wing allows a reasonable HI fit that 
classifies this absorber as a sub-DLA. 

For the z abs = 2.9422 sub-DLA towards QSO A we de- 
tect a range of metal species: Si II, Fe II, Zn II, S II, Ni II, 
Al II and Al III, see Table [5] and Figure [7] N V A 1242 may 
be present, but is weak (peak optical depth r ~ 0.1) and 
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Table 3. Ionic Column Densities for £ a b s — 2.6617 DLA towards QSO A 



Ion 


A 


/-value 


log N 


log N adopt 


H I 


1215.6701 


0.41640 




20.48±0.10 


Al III 


1854.7164 


0.5390 


12.97±0.10 


12.97±0.10 


Al III 


1862.7895 


0.2680 


12.99±0.14 




Fe II 


1608.4511 


0.0580 


14.41±0.05 


14.36±0.10 


Fe II 


1611.2005 


0.001360 


<14.5 




Fe II 


2344.2140 


0.1140 


14.22±0.03 




Fe II 


2374.4612 


0.0313 


14.30±0.06 




Fe II 


2382.7650 


0.3200 


14.12±0.03 




Si II 


1808.0130 


0.002186 


15.05±0.11 


15.05±0.11 


Zn II 


2026.1360 


0.4890 


12.60±0.13 


12.40±0.20 


Zn II 


2062.6640 


0.2560 


12.40±0.20 




Cr II 


2056.2539 


0.1050 


<12.7 


<12.7 


Ni II 


1741.5531 


0.04270 


13.35±0.20 


13.35±0.20 



All limits arc 3a. 



Table 4. Ionic Column Densities for 2 a b s — 

2.6625 sub-DLA towards QSO B 



Ion 


A 


/-value 


log N 


log N adopi 


H I 


1215.6701 


0.41640 




20.11±0.10 


Fc II 


2344.2140 


0.1140 


13.54±0.09 


13.87±0.13 


Fe II 


2374.4612 


0.0313 


13.87±0.13 




Si II 


1808.0130 


0.002186 


<14.4 


<14.4 


Zn II 


2026.1360 


0.4890 


<12.1 


<12.1 


Cr II 


2056.2539 


0.1050 


<12.7 


<12.7 


Ni II 


1741.5531 


0.04270 


<13.0 


< 13.0 



All limits arc 3<r. 



Table 5. Ionic Column Densities for z abs = 2.9422 sub-DLA towards QSO A 



Ion 


A 


/-value 




log N 




log N adopt 


H I 


1215.6701 


0.41640 








20.28±0.05 


S II 


1250.5840 


0.005453 


15 


,01±0. 


10 


15.01±0.10 


S II 


1253.8110 


0.01088 


15 


,01±0. 


08 




Al II 


1670.7874 


1.880 


13 


.84±0. 


06 


^13.84 


Al III 


1854.7164 


0.5390 


13 


73±0. 


04 


^13.73 


Al III 


1862.7895 


0.2680 


13 


.73±0. 


07 




Fe II 


1608.4511 


0.0580 


14 


69±0. 


04 


14.69±0.04 


Fe II 


1611.2005 


0.001360 


14 


66±0. 


21 




Si II 


1808.0130 


0.002186 


15 


,34±0. 


08 


15.34±0.08 


Zn II 


2026.1360 


0.4890 


12 


.62±0. 


12 


12.40±0.33 


Ni II 


1370.1310 


0.07690 


13 


,71±0. 


23 


13.66±0.28 


Ni II 


1454.842 


0.0323 


13 


60±0. 


21 





All limits are 3a. 



Table 6. Ionic Column Densities for 2 a bs — 2.9467 sub-DLA towards QSO B 



Ion 


A 


/-value 


log N 


log N adopt 


H I 


1215.6701 


0.41640 




19.70±0.10 


S II 


1253.8110 


0.01088 


<14.0 


<14.0 


Al II 


1670.7874 


1.880 


12.57±0.10 


12.57±0.10 


Al III 


1862.7895 


0.2680 


12.51±0.22 


< 12.6 


Fe II 


1608.4511 


0.0580 


13.29±0.25 


13.29±0.25 


Si II 


1808.0130 


0.002186 


<14.6 


<14.6 


Zn II 


2026.1360 


0.4890 


<12.1 


<12.1 


Ni II 


1370.1310 


0.07690 


<13.4 


<13.4 



All limits arc 3<r. 
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blended with another feature. We do not quote a detection 
or limit for Cr II since all the lines are in a contaminated part 
of the spectrum. Al II A 1670 is clearly saturated, so we only 
quote a lower limit for N(A1 II). For Al III AA 1854, 1862 
the two lines give consistent column densities, despite their 
high equivalent widths. However, to be cautious, we quote 
their column densities as lower limits. The Fe II A 1608 line 
towards QSO A is quite strong and may be slightly satu- 
rated. However, the much weaker Fe II A 1611 is marginally 
detected and gives a consistent column density. 

We have the same problem for Zn II as for the DLA 
at Zabs ~ 2.66 discussed in the previous subsection, i.e. the 
potential contribution from Mg I to the Zn II A 2026 line. 
However, the situation is even more complicated in this case 
because not only is Mg I A 2852 not covered, but also the 
Zn II A 2062 line is in a highly contaminated part of the 
spectrum, so can not be used and we have no estimate of the 
Cr II contribution. We therefore estimate an upper limit to 
the contribution of Mg I by considering values measured in 
other metal-rich absorbers. From Peroux et al. (2006) and 
Herbert-Fort et al. (2006), we assume a rest-frame upper 
limit contribution of 50 mA from Mg I. This means that 
N(Zn II) may require a downward revision of up to 0.4 dex. 
We therefore quote a column density of Zn II with error 
bars that account for this range of possibilities, see Table 
[5] Regardless of our uncertainties in N(Zn), this absorber 
is clearly metal rich, since the column densities of Si and S 
also yield abundances ~ 1/3 Zq. 

We also search for H2 in the sub-DLA towards QSO A. 
At z a bs = 2.9422 we have coverage of Lyman J = 0, 1 rota- 
tional bands for 6 — to — vibrational transitions. These 
first two J states usually dominate the molecular column 
density in DLAs (e.g. Ledoux et al. 2003). Assuming a red- 
shift matched to the central position of the metal lines, i.e. 
z a bs = 2.9422, we combine the limits on J — and J = 1 to 
determine a limit of log N(Ha) < 14.5 in this line of sight. 

For the sub-DLA towards QSO B, we only determine 
column densities for Fe II and Al II and upper limits for the 
other species listed in Table [(5] Al III A 1854 is in a region 
of moderate contamination in the B spectrum, so that the 
column of a weak line is unreliable. There is a marginally 
significant (2.6<r) feature close to the expected position of 
Al III A 1862, but since it is below the 3a level, we quote it 
as a limit. 



3 DISCUSSION 
3.1 Abundances 

The final abundances for the three absorbers with log N(HI) 
> 20 are given in Table [7] The 2 absorbers in QSO A both 
have very high metallicities, based on the abundance of Zn: 
1/3 and 1/5 Zq. Such high metallicities are very rare even in 
low redshift DLAs (e.g. Meiring et al. 2006), leading to the 
widespread conclusion that the DLA cross-section is dom- 
inated by metal-poor gas. It has been suggested that this 
may be due to a dust-induced bias against metal-rich galax- 
ies, although there is currently no observational evidence 
to support this claim (e.g. Ellison et al. 2001, 2004, 2005; 
Akerman et al. 2005; Jorgenson et al. 2006). Herbert- Fort 



et al. (2006) have shown that selecting absorbers from the 
SDSS on the basis of strong metal lines can readily identify 
DLAs with metallicities in excess of O.IZq, indicating that 
metal-rich systems do exist at high redshift, but are simply 
rare. The two DLAs studied here have metallicities consis- 
tent with the metal-strong population of Herbert-Fort et al. 
Moreover, the abundances that we determine are compara- 
ble to those measured in cB58 (2/5 Zq at z = 2.7, Pettini et 
al. 2002) and other LBGs at 2 < z < 3 (Teplitz et al. 2000; 
Shapley et al. 2004). 

We also determine relatively high depletion factors for 
the DLAs based on [Zn/Fe] = +0.88, +0.55 in the DLAs 
at Zabs = 2.66 and 2.94 respectively. These values are usu- 
ally considered as measures of the dust:metals ratio, since 
these two elements trace each other well in most Galactic 
stars (although see the caveats in Nissen et al. 2004), yet 
Fe is highly refractory whereas Zn is not (Savage & Sem- 
bach 1996). The relative abundances measured here imply 
that ~ 15 - 30 % of the metals in these two DLAs are in the 
gas phase and follows the broad trend of increased depletion 
with increasing metallicity (e.g. Meiring et al. 2006 for the 
most recent compilation). Although quite extreme by DLA 
standards, these depletion factors still barely overlap with 
measurements in the Galactic disk and Magellanic Clouds 
(e.g. Roth & Blades 1997). 

Although this sightline intersects gass associated with 
the diffuse ISM, given its high metallicity (1/3 Zq) and high 
depletion, it is perhaps surprising that we do not detect 
H2 down to a molecular fraction of log /(H2) = 2N(H2) / 
[2N(H 2 ) + N(HI)] < -5.5 for the z abs = 2.94 DLA0. Ledoux 
et al. (2003) and Petitjean et al. (2006) have suggested that 
both metallicity and depletion may be important factors 
for a galaxy to be able to form (and maintain) a signifi- 
cant column density of molecules. As a comparison with the 
Zabs = 2.94 DLA studied here, we can consider molecular 
fractions in LMC sightlines, where the metallicity is similar 
to this DLA. Tumlinson et al. (2002) showed that H2 is de- 
tected in ~ 50% of LMC sightlines, so we may naively expect 
a similar detection rate in other galaxies with similar metal- 
licities. Although the statistics for metal-rich DLAs are still 
poor, a 50% detection rate does broadly fit the high redshift 
data for DLAs with metallicities above 1/5 Zq (Petitjean et 
al. 2006). However, metallicity is clearly not the only factor, 
since the SMC, whose metallicity is a factor three lower than 
the LMC, has detections for 90% of the sightlines studied by 
Tumlinson et al. (2002) . This has been proposed to be due to 
higher star formation rates in the LMC which preferentially 
photo-dissociate molecules. 

Finally, we note that the z a bs = 2.94 DLA has a relative 
abundance of S and Zn that is close to the solar value, [S /Zn] 
= +0.05. In Galactic stars, we observe enhanced a/Fe ratios 
at metallicities below [Fe/H] < —1.0, due to the time delay 
between SNII and SNIa enrichment. In DLAs, the two unde- 
pleted elements S (an a element) and Zn (which traces the 
Fe-peak in Galactic halo stars at [Fe/H] > —2.5) are usu- 
ally preferred tracers of the relative contributions of SNIa 



3 Although it is possible that we have under-estimated the con- 
tribution of Mg I to the Zn II A 2026 line and the true metallicity 
may be less than we have deduced. However, the abundance based 
on Si or S is also X/3Zq. 
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Table 7. Abundances for Absorbers with N(HI) > 20.1 



QSO 


^abs 


Log N(HI) 


[Fe/H] 


[Si/H] 


[Zn/H] 


[Cr/H] 


[S/H] 


[Ni/H] 


A 


2.6617 


20.48±0.10 


-1.59 ±0.14 


-0.97 ±0.15 


-0.71 ± 0.22 


< -1.43 




-1.35 ± 0.22 


B 


2.6625 


20.11±0.10 


-1.71 ±0.16 


< -1.25 


< -0.51 


< -1.06 




< -1.33 


A 


2.9422 


20.28±0.05 


-1.06 ±0.06 


-0.48 ±0.09 


-0.51 ± 0.33 




-0.46 ± 0.11 


-0.84 ±0.28 



All limits arc 3o\ 



and SNII (e.g. Nissen et al. 2004). The universally low (typi- 
cally sub-solar) values of [S/Zn] in DLAs indicates that their 
chemical enrichment history is very different to that of the 
Milky Way. New measurements from local dwarf galaxies in- 
dicate that low a/Fe ratios may actually be common in the 
Local Group, e.g. Shetrone et al. (2003). At a metallicity of 
1/3Zq the value of [S/Zn] measured here is actually inter- 
mediate between typical Galactic values (e.g. Gratton et al. 
2003) and those found in the more metal rich Local Group 
dwarfs such as Fornax, the LMC and Sagittarius (Pompeia, 
Hill & Spite 2005; Bonifacio et al. 2004). 



3.2 Are DLA galaxies over 100 h 7( j kpc in size? 

At low redshift, we have a number of measurements of DLA 
absorber size. At z = 0, Zwaan et al. (2005) showed that 
a 21cm survey of selected UGC galaxies finds typical di- 
mensions for DLA-equivalent column densities ranging from 
~ 5 - 50 kpc, depending on galaxy mass and orientation. 
At similarly low redshifts (cz < 4000 km/s) Bowen, Pettini 
& Blades (2002) found that if a line of sight passes within 
200 kpc of a nearby galaxy, the covering factor for HI gas 
with log N(HI) > 13 is effectively 100%. However, the cross- 
section of gas that is optically thick at the Lyman limit will 
be considerably smaller than for absorbers with 13 < log 
N(HI) < 16. At these lower column densities, the HI absorp- 
tion can presumably be associated not only with extended, 
diffuse halos, but also with non-galactic structures such as 
filaments and intra-group media. At intermediate redshifts, 
Chen & Lanzetta (2003) determined a characteristic size 
for DLAs of ~ 25 kpc, and also found that a significant 
number of DLA fields have multiple galaxies at the absorber 
redshift. In contrast, constraints on galaxy/absorber size on 
kpc scales at z > 2 are very few. They are limited to 2 pairs 
of QSO sightlines which probe transverse scales of 5-10 kpc 
(Lopez et al. 2005; Smette et al. 1995), a measurement of 
spatially offset Lya fluorescence (Adelberger et al. 2006), a 
handful of extended radio sources (Foltz et al. 1988; Briggs 
et al. 1989) and extended [OIII] emission from 3 DLAs at 
2 < z < 3 (M0ller, Fynbo & Fall 2004; Weatherley et al. 
2005). Although these observations generally imply gaseous 
extents of at most a few tens of kpc, Labbe et al. (2003) 
have used deep imaging to detect a population of galaxies 
at 1.5 < z < 3 whose disks have a comparable size to the 
Milky Way and Prochaska & Wolfe (1997) have also inter- 
preted DLAs as large disks. The critical question is there- 
fore: is it likely that the common absorption seen in SDSS 
1116±4118 AB is due to galaxies with 100 kpc gas halos at 
z ~ 3? 

Some first clues as to whether the coincident absorp- 
tion is due to a common structure can be gleaned from 
the relative velocities across the line of sight. Based on the 



metal lines, the velocity separations between the DLA/sub- 
DLA in QSO A/B are 65 km/s (z abs ~ 2.66) and 340 km/s 
(«abs ~ 2.94). Although only a handful of rotation curves 
and velocity fields have been studied at high redshift (e.g. 
Erb et al. 2003, 2004), a significant number have velocity dis- 
persions in excess of 200 km/s and rotation velocities > 150 
km/s (i.e. a total velocity difference > 300 km/s). There- 
fore, although the shear between QSO A/B is large for the 
z a bs ~ 2.94 absorber, it is still feasible that it is within the 
same galaxy. 

We can also appeal to the metallicities for clues to the 
absorbers' structure. At both z ~ 2.66 and z ~ 2.94, at least 
one of the absorbers has a high metallicity. On the one hand, 
the existence of a mass-metallicity relation at all redshifts 
studied (Tremonti et al. 2004; Savaglio et al. 2006; Erb et al. 
2006) indicates that these galaxies may have stellar masses 
at least ~ 10 9 — 10 9 ' 5 Mq, possibly lending support to the 
idea of a large physical extent. On the other hand, such high 
abundances seen in absorption, consistent with metallicities 
inferred from the emission line measurements in LBGs, in- 
dicates that these sightlines may intersect the central part 
of their galaxies. If true, it would imply even larger sizes 
for the absorbing galaxy than the simple interpretation of 
common coincidence, i.e. r > 100 h^Q kpc, rather than d > 
100 ZiyQ 1 kpc. However, neither the information on metal- 
licity or velocity provide very satisfactory conclusions as to 
the nature of the coincident absorption. We therefore turn 
to galaxy simulations to help us interpret our observations. 



3.2.1 Cosmological simulations of galaxy formation 

The code used for the simulations is a significantly improved 
version of the TreeSPH code, which has been used previously 
for galaxy formation simulations (Sommer-Larsen, Gotz & 
Portinari 2003). The main improvements over the previous 
version are: (1) The 'conservative' entropy equation solving 
scheme (Springel & Hernquist 2002) . (2) Non-instantaneous 
gas recycling and chemical evolution, tracing 10 elements 
(Lia et al. 2002a,b); the algorithm includes supernovae of 
type II and type la, and mass loss from stars of all masses. 
(3) Atomic radiative cooling depending both on the metal 
abundance of the gas and on the meta-galactic UV field, 
modeled after Haardt & Madau (1996) is invoked, as well as 
simplified treatment of radiative transfer, switching off the 
UV field where the gas becomes optically thick to Lyman 
limit photons on scales of ~ 1 kpc. 

In this paper, very high resolution simulations of two 
galaxies, known to become large disk galaxies at 2=0, are 
analysed in relation to DLA and sub-DLA sight-line prop- 
erties. We take the two galaxies simulated by Razoumov & 
Sommer-Larsen (2006), K15 and K33, which have, at 2=0, 
characteristic circular speeds of V c =245 and 180 km/s, re- 
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Figure 8. Projected N(HI) map for the central 200 kpc for our two simulated galaxies, K15 (left) and K33(right) at Z = 3.6,3.0 and 
2.3 (top through bottom panels). This projection is through the XY (face-on) plane. The bar in the top right panel shows the typical 
transverse separation at the absorber radshift for the QSO pair studied here. The white bar in the top right panel shows 110 fhfg kpc, 
the approximate physical separation of the 2 lines of sight at the absorbers' redshift in SDSS1116+4118 AB. 
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Figure 9. Projected N(HI) map for the central 60 kpc for our two simulated galaxies, K15 (left) and K33(right) at z = 3.6,3.0 and 2.3 
(top through bottom panels). This projection is through the XY (face-on) plane. 
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Figure 10. Projected N(HI) map for the central 60 kpc for our two simulated galaxies, K15 (left) and K33(right) at z = 3.6, 3.0 and 2.3 
(top through bottom panels). This projection is through the XZ (edge-on) plane. 
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K33 z=0.3 XY 




Figure 11. Projected N(HI) map for the lowest redshift slice currently simulated for galaxy K33. Two scales and two orientations are 
shown, as for Figures 151 to ITul The specific angular momentum, j, of the galaxy is 660 km s _1 kpc, near the median of the observed 
range of j ~ 350 — 1400 km s — 1 kpc typical for a disk galaxy of V c =180 km s — 1 at z ~ (e.g., Sommer-Larsen et al. 2003). 



spectively. The galaxies bracket typical disk galaxy forma- 
tion histories: the formation of the larger K15 disk is merger 
induced (e.g., also Robertson et al. 2004), with the disk 
growing strongly between 2=1 and 0, whereas K33 already 
starts developing a disk by z ~ 2.5, which subsequently 
grows gradually to the present epoch. The z=Q virial masses 
of the two galaxies are M vir = 8.9xlO n and 3.7xW n M Q , for 
K15 and K33 respectively. 

The two galaxies were simulated in the standard ACDM 
cosmology using the 'zoom-in' technique (e.g., Sommer- 
Larsen 2006) to study the formation and evolution of in- 
dividual galaxies in full cosmological context. Total particle 
numbers used were 2.2 and 1.2 million, for K15 and K33 
respectively. Particle masses and gravitational (spline) soft- 
ening lengths of m gas =m»=8.45xl0 4 and rriDM=4.83xl0 5 



h^Mo, and e gas =e,=176 and eoM=316 h^pc, respec- 
tively, were adopted. The gravity softening lengths were 
fixed in physical coordinates from z~6 to z~0, and in co- 
moving coordinates at earlier times. The minimum SPH 
smoothing length in the simulation was about 12/ifg 1 pc. A 
Kroupa IMF (Kroupa 1998) was used in the simulations, and 
early rapid and self-propagating star-formation (sometimes 
dubbed 'positive feedback') was invoked (Sommer-Larsen et 
al. 2003). 

Although K15 and K33 galaxies have been chosen from 
our simulation to represent two different disk galaxy evo- 
lutionary paths (and end products) there are some generic 
features to both galaxies. For example, they both tend to be 
larger and have a more extended gas distribution at higher 
redshift. At low redshift, the gas has cooled and settled into 
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a more organised and centrally concentrated structure. How- 
ever, proto-disks can form even at quite high redshifts (see 
the top left panel of Figures [8] and [9}. The two different 
scales for these figures highlight that the majority of DLA 
gas is distributed on scales of tens, rather than hundreds, 
of kpc. Moreover, the region of the galaxy which is actively 
star formating is very small (~ kpc scale) compared with 
the distribution of high column density gas. This is quali- 
tatively consistent with the observation of low abundances 
in DLAs (e.g. Pettini et al. 1999) compared with emission 
line abundances in absorption selected galaxies (e.g. Ellison, 
Kewley & Mallen-Ornelas 2005; Chen, Kennicutt & Rauch 
2005). The spatial separation of the bulk of DLA-inducing 
gas from the regions of high star formation has also been 
inferred by the lack of low surface brightness galaxies in the 
Hubble ultra deep field (Wolfe & Chen 2006). 

We analyse the model galaxies by passing sightlines 
through the simulation box and integrating the HI volume 
density along the line of sight. In Figures [51 to 1101 we show 
two projections of the two simulated galaxies where each 
box has been projected into a 600 x 600 grid. The colour 
contours show the projected HI gas column density in units 
of atoms cm -2 and are therefore analogous to the observed 
N(HI) columns measured in QSO sightlines. The N(HI) at 
each point along the line of sight is evaluated by weight- 
ing all of the particles within the local smoothing length h 
(where 2h is the distance within which 50 particles are lo- 
cated) according to the spline smoothing kernel (Monaghan 
& Lattanzio 1985). For this work, we use two different ren- 
derings of each galaxy in order to alleviate orientation ef- 
fects. These two projections correspond to projecting lines 
of sight through the XY and XZ faces of the same simula- 
tion cube. The XY plane corresponds roughly to a face-on 
configuration of the central star-forming disk whereas the 
XZ is approximately edge-on. However, as can be seen from 
a visual inspection of Figures ISlto 1101 the distribution of HI 
gas is not significantly different even in these two 'extreme' 
orientations; this is checked quantitatively below. Although 
we do not use simulations at lower redshifts, for demonstra- 
tion pruposes we include here images of the K33 at z = 0.3, 
the lowest redshift to which the simulation has currently 
progressed. Figure [TT] shows the galaxy on two spatial scales 
and in both the face-on and edge-on orientations. These im- 
ages show the stability and growth of the disk since z = 3.6. 
The K15 simulation has not yet progressed beyond z = 2.0, 
so we can not show a lower redshift rendition. 

In order to calculate the probability of finding a DLA 
or sub-DLA in each of a pair of sightlines, we construct a 
sample of 200,000 pairs from each of the simulated galaxies 
at 2 = 2.3, 3.0 and 3.6. This suite includes 100,000 sightlines 
from each of the face-on and edge-on renditions. In every 
pair, the first sightline is drawn at random from the sub- 
set known to intersect sub-DLA (log N(HI) > 19.5) or DLA 
(log N(HI) 20.3) gas. The second sightline is drawn at 
random from the full simulation. In Figure Q2] and [T3] we 
show the fraction of pairs with a sub-DLA or DLA in the 
second line of sight, given the presence of a sub-DLA or DLA 
in the first for galaxies K15 and K3J0. For example, for the 
larger K15 galaxy at z = 3.0 about 6% of sightlines with a 



4 In this simulation we classify sub-DLA as absorbers with log 



DLA in the first line of sight also have a DLA in the second 
for separations 20 - 30 kpc. For sub-DLAs, this fraction 
increases to ~ 10%. At separations from 40 to 150 kpc, the 
likelihood of finding coincident DLA absorption decreases 
from ~ 1% to 0.1% in K15. In the smaller K33 galaxy, these 
probabilities are about a factor of 5 smaller at z = 3.0 and 
z = 3.6. 

In our observations of SDSS 1116+4118 AB we have 
one case of a sub-DLA with no counterpart in the second 
line of sight (z a b s = 2.47), one case of a DLA with sub- 
DLA counterpart (z a h s = 2.66) and one case of a pair of 
sub-DLAs (zabs = 2.94) all with separations ~ 110 /ifg 1 kpc. 
From Figures [T2l and H"3"l it can be seen that the probability 
of these two 'double hits' is < 10 -5 . The present simulations 
still fail to reproduce the observed angular momenta of ac- 
tual disk galaxies by about a factor of 1.5, and hence the real 
(proto-)disk galaxies may be up to 50% more extended than 
the present models. However, the probability of a double in- 
tersection is still < 1% at 110/1.5 ~ 75 h^g kpc. Another 
potential uncertainty in the models is our simplified treat- 
ment of radiative transfer, which governs the transition from 
neutral to ionized gas and can therefore affect the effective 
area of projected sub-DLA absorbers. Experiments with our 
radiative transfer treatment indicate that this effect is likely 
to affect sub-DLA numbers by, at most, a factor of two. 
Simulations that contain a more complete radiative transfer 
calculation are currently underway and will be presented in 
a future paper (see also Razoumov et al. 2006). Even with 
the uncertainties of angular momentum and the transition 
between neutral and ionized gas, the conclusion that a single 
large galaxy, such as those tested here, is unlikely to cause 
coincident absorption in SDSS 1116+4116 is unchanged, al- 
though tests with a larger suite of galaxies would be desir- 
able. As the two simulation volumes also include satellites 
and merging components around the (proto-)galaxies (e.g. 
bottom left panel of Figure [SJ , the low probability of coin- 
cident pairs found above also argues against a high covering 
fraction of satellite systems or an actively merging galaxy 
(such as that observed by Miley et al. 2006) as the explana- 
tion for the common absorption across the line of sight, at 
least based on these two simulated galaxies. Again, a wider 
suite of simulations is needed to confirm this for a more 
extensive parameter space. 



3.2.2 DLA clustering analysis 

Having rejected a large single system, and a galaxy plus 
satellite structure as the reason for coincident absorption 
towards SDSS 1116+4118 AB, we now consider the possibil- 
ity that we are observing a structure that contains multiple 
galaxies. Such a possibility is perhaps not completely un- 
expected given that galaxy groups represent the most com- 
mon environment at least at low redshift (e.g Tully 1987). 
Galaxy clustering at z ~ 3 is becoming increasingly well- 
documented, such as the extensive literature on LBGs e.g. 
Giavalisco et al. (1998), Steidel et al. (1998), Porciani & 
Giavalisco (2002), Adelberger et al. (2005a,c), Steidel et al. 



N(HI) 19.5 gas, i.e. grouping together classical DLAs and sub- 
DLAs into the same category. 
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Figure 12. The fraction of second sightlines in a pair, as a function of separation, that will exhibit a DLA (red) or sub-DLA (black) 
given the presence of a DLA or sub-DLA in the first. This probability function for finding common absorption is calcuated for 3 different 
redshifts (see panel labels) for the more massive galaxy in our simulation, K15. The top three panels are for the XY (face-on) orientation, 
and the bottom three panels are for the XZ (edge-on) orientation. The pair probabilities for the two orientations are very similar, as may 
be expected from a visual inspection of Figures 151 to 1101 
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Figure 13. As for Figure IT2l but for the less massive galaxy in our simulation, K33. 
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Figure 14. Probability of coincident absorption as a function of 
column density threshold. The lower dotted line shows the cos- 
mic average probability for finding an absorber within a ±400 
km/s window at redshift z = 2.7 in the absence of clustering. 
The solid line represents the probability for finding an absorber 
in sightline B, 13.8" away from and within ±400 km/s of an ab- 
sorber detected in sightline A (or vice versa). This prediction 
assumes that high-column density absorbers have a power law 
auto-correlation function £ = (r/fo)'— 7) which is the same as the 
DLA-LGB correlation function measured by Cooke et al. (2006) 
(7 = 1.6; ro = 2.66^'g), and the blue-shaded region represents 
their ±1<t errors. The dashed illustrates the effect of a steeper 
correlation function (7 = 2.1; ro = 2.8) which was found to be 
consistent with a subset of the Cooke et al. (2006) data. 



(2005). Although these investigations have selected galax- 
ies based on their emission properties, Cooke et al. (2006) 
have shown that DLAs have similar spatial distributions as 
the LBGs. DLAs may therefore also be highly clustered at 
z ~ 3. This posit is supported by recent observations of 
Zibetti et al. (2007) who stack images of the QSO fields of 
SDSS-selected Mg II absorbers. They find extended starlight 
out to ~ 200 kpc, indicating that there can be galaxies at the 
redshift of the absorber out to these impact parameters, con- 
sistent with a group environment. Chen & Lanzetta (2003) 
have also found that the fields of DLAs at z < 1 often have 
multiple galaxies at the absorber redshift. We can estimate 
the likelihood of a multiple galaxy coincidence by consider- 
ing the clustering scale of DLA galaxies and their column 
density distribution. 

In the absence of clustering, the line density of ab- 
sorption line systems per unit redshift above the column 
density threshold N(HI) is given by the cosmic average 
(^gr) (>N(HI),z). At an average location in the universe, the 
probability of finding an absorber in a background quasar 
spectrum within the redshift interval Az = 2(1 + z)Av/c, 
corresponding to a velocity interval ±Av is simply P = (dN / 
dz)Az. For a close pair of quasar sightlines, the presence of 
an absorber in one of the sightlines implies an increased 
probability of finding an absorber at a similar redshift in 
the neighbouring sightline. If the quasar sightlines have a 
transverse separation R, and assuming that one searches a 
velocity interval ±Av about the absorber redshift, we follow 



Hennawi & Prochaska (2007) and express the enhanced line 
density in the neighbouring sightline, ^-(R, Av), in terms 
of a transverse correlation function \± (R) as 

dN ,„ „ . /dN\ , ,„ A 

where x±(R>Av) is given by an average of the 3-d 
absorber auto-correlation function, f(r), over a cylindrical 
volume with cross-section A equal to the absorption cross 
section, and height given by the length in the line-of-sight 
direction corresponding to the velocity interval 2Av. Pro- 
vided that we are in the 'far-field' limit, where the transverse 
separation is much larger than the absorber cross section 
R 3> y/A, it is a very good approximation to replace the 
volume integral over the cylinder with a line integral over 
the range L = 2Av / aH (z), where a is the scale factor and 
H(z) is the Hubble constant at redshift z (see Hennawi & 
Prochaska 2007 for details). 

Studies of the clustering of high-column density ab- 
sorbers, such as DLAs, are hampered by the relatively small 
number of absorber pairs known, which is in turn limited by 
the number of quasar pairs with separations smaller than 
the absorber auto-correlation length. These poor statistics 
can be circumvented if instead one cross-correlates absorbers 
with galaxies. Cooke et al. (2006) recently measured the 
clustering of LBGs around DLAs at z ~ 3 and measured a 
best fit cross-correlation length of ro = 2.66^'g h" 1 Mpc 
assuming a fixed power-law slope of 7 = 1.6 (see also Adel- 
berger et al. 2003; Gawiser et al. 2001; Bouche & Lowenthal 
2004) . This clustering measurement is very close to the auto- 
correlation of LBGs, indicating that the auto-correlation of 
DLAs is similar in strength. 

In Figure Q3] we use the Cooke et al. (2006) measure- 
ment and equation [T] to estimate the probability of detect- 
ing coincident high-column density absorbers as a function 
of column density threshold. The lower dotted line shows 
the cosmic average probability P = {dN / dz) Az for finding 
a single absorber within a ±400 km/s window at redshift 
2 = 2.7 in the absence of clustering. Calculating this prob- 
ability requires an expression for the column density distri- 
bution f(N), which we obtained by joining the the f(N) for 
DLAs log N(HI) ^ 20.3 (Prochaska et al. 2005) with that 
for sub-DLAs 19 < log N(HI)< 20.3 (O'Meara et al. 2007). 
The solid line represents the probability for finding an ab- 
sorber above a given N(HI) in a second sightline, given the 
presence of an absorber in the first sightline within a ±400 
km/s velocity window. This calculation uses the Cooke et 
al. (2006) best-fit measurement for 7 = 1.6 and is performed 
for the specific case of a 13.8 arcsec separation. The blue- 
shaded region represents the range permitted by the ±lcr 
errors in 7 quoted by Cooke et al. (2006). For this model, 
the probability of finding coincident sub-DLA absorption 
(both absorbers N(HI) > 19.5) in the SDSS 1116+4118 AB 
binary within a ±400 km/s window is about 3%. In the ab- 
sence of clustering, this value is about 0.5%, i.e. the random 
probability of finding a sub-DLA in a ±400 km s" 1 velocity 
window based on the number density distribution alone. 

The smallest separation probed by the Cooke et al. 
(2006) measurement is (co-moving) ~ 400 h^ Q kpc, or about 
40% larger than the separation of our quasar pair sightline. 
The auto-correlation functions of high redshift galaxies tend 
to become progressively steeper on (co-moving) scales less 
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than about 1 h^ 1 Mpc, characteristic of the sizes of dark 
matter halos (Coil et al. 2006; Ouchi et al. 2005; Lee et al. 
2006; Conroy et al. 2005), and similar behaviour might be 
expected around the absorbers considered here. To illustrate 
the effect of a steeper correlation function, the dashed line 
in Figure [14] shows the probability for coincident absorp- 
tion for a correlation function with 7 = 2.11 and ro = 2.81. 
This choice was motivated by the fact that Cooke et al. 
(2006) found that a subset of their data (11 of the 15 DLA 
sightlines) were best-fit with cross-correlation function pa- 
rameters 7 = 2.11^j'4 and ro = 2.8l2.o- F° r this model, the 
probability of finding coincident sub-DLA absorption (both 
absorbers N(HI) > 19.5) in the SDSS 1116+4118 AB binary 
within a 400 km/s window is about 8%. 

Although the probabilities that we determine from the 
clustering analysis are not large (< 10% for a single coinci- 
dent pair), they are more than an order of magnitude more 
likely than the probability of absorption from a single large 
galaxy (previous sub-section). Previous observations of some 
wide separation QSO groupings have also concluded that in- 
tervening (super-)clusters may be responsible for correlated 
absorption (e.g. Jakobsen et al. 1986, 1988; Francis & Hewett 
1993), although the scales probed here are an order of mag- 
nitude smaller than those studies. We note that one of the 
coincident absorbers is separated from the systemic redshift 
of the QSO by ~ 3000 and 4500 km/s (z abs = 2.94 in QSO 
A/B respectively). Such 'proximate' DLAs (PDLAs) have 
been shown to exhibit excess clustering around QSOs (El- 
lison et al. 2002; Russell, Ellison & Benn 2006; Prochaska, 
Hennawi & Herbert-Fort 2007). The probability for coinci- 
dent absorption close to the QSO redshift may therefore be 
further enhanced beyond the calculation above. 

Finally, we note that these clustering statistics are not 
affected by any a priori knowledge of the presence of a DLA 
in either sightline. The analysis assumes the presence of one 
DLA and calculates the probability that a matching DLA 
is found in the second sightline. Of course, all pairs with at 
least one DLA (i.e. regardless of whether there is a coinci- 
dent absorber in the second sightline) should be included 
in order to calculate unbiased statistics. Currently, SDSS 
1116+4118 AB is our only fully analysed example, although 
we discuss in the final section the prospects for improving 
these statistics. 



3.3 Implications for line of sight observations 

Regardless of the physical interpretation of the coincident 
absorption (single large galaxy, clustering etc.), the corre- 
lated absorption has interesting ramifications for single line 
of sight observations. The most simplistic extrapolation of 
the three intervening absorbers towards this one binary sys- 
tem indicates that ~ 2/3 of DLAs are not simple systems. 
From a different spatial perspective, these spatially resolved 
systems could be superimposed along a single sightline lead- 
ing to an absorber that is a composite of two (or more) 
components. Whilst the consideration and treatment of com- 
plex absorption systems consisting of many super-imposed 
'clouds' is not new (e.g. Jenkins 1986; M0ller, Jakobsen & 
Perryman 1984), the data presented here supply the first ob- 
servational suggestion that superpositions of multiple galax- 
ies may be relevent. For the velocity differences observed 
here (65 - 340 km/s) the Lya would appear as a single 
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Figure 15. Illustration of the potential effect of line of sight 
blending; Lya absorption is shown in the left-hand panels and 
Zn II A 2026 is shown in the right-hand panels. In all cases, the 
labels within the panels show the input parameters of the two 
individual simulated absorbers and the labels above the panels 
indicates the value that would be measured for the combination 
of the two. 



trough with a column density indistinguishable from the sum 
of the components. Indeed, only DLAs whose relative veloc- 
ities are above ~ 1000 km/s would be recognised as multiple 
absorbers from the Lya line alone (depending on the column 
densities, since higher N(HI) systems need to be separated 
by larger velocities for identification) . Although a two DLA 
blend could theoretically be identified by the asymmetry 
between the blue and red wings of the profile, in practice 
the delicate process of continuum fitting could obscure this. 
Difficulties include the complex underlying form of the QSO 
continuum (such as emission lines), flux suppression in the 
Lya forest and, in the case of cross-dispersed spectra, com- 
plexity of removing the blaze, joining the orders and poor 
flux calibration. Indeed, when fitting the continuum around 
the Lya trough of a DLA, curvature is often required to yield 
a normalized spectra that can be well-fitted with a damped 
profile (e.g. Prochaska et al. 2003) , which could lead to over- 
sight of a blended absorber. 

Cases of (partially) resolved multiple DLAs have pre- 
viously been reported by several authors (e.g. Wolfe et 
al. 1986; Ellison & Lopez 2001; Lopez & Ellison 2003; 
Prochaska et al. 2003). Peroux et al. (2003) have also noted 
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that 30% of their sub-DLA sample are associated with an- 
other absorber. The impact of superimposed absorbers de- 
pends on a variety of factors, such as relative velocities, 
N(HI) and metallicities. We illustrate three possible combi- 
nations in Figure [Pol In the first scenario, a DLA is blended 
with a sub-DLA, both of which have the same metallicity, 
but due to the lower N(HI) of the sub-DLA, its metals are 
below the detection threshold of the spectrum, leading to 
an under-estimate of the abundance in the blended sightline 
(top panel, Figure [To)) . A similar under-estimate would re- 
sult if the velocity difference is large enough that the metal 
lines from the second system are not included in the fit 
(middle panel, Figure I15[) . In practice, however, this par- 
ticular situation is unlikely to occur since normally other, 
stronger, metal lines such as Fe II would signal the existence 
of high velocity gas. Alternatively, if a blended sub-DLA had 
a higher intrinsic abundance, as has been proposed to be 
systematically the case (e.g. Peroux et al. 2006) we would 
infer an artificially high metallicity for the DLA (bottom 
panel, Figure I15[) . Dessauges-Zavadsky et al. (2006) have 
noted that the abundance ratios of different velocity com- 
ponents in single line of sight DLAs are often non-uniform, 
which is expected if the system is a blend of more than one 
galaxy. These effects are of course most severe when the 
N(HI) values of the blended absorbers are close in value; a 
sub-DLA with log N(HI) = 19.0 will have little impact on a 
DLA with log N(HI) = 21.0. 

The superposition of DLAs in single lines of sight may 
partially explain the often complex nature observed in metal 
line profiles and the extreme rarity of DLAs with simple ve- 
locity structures. The potential for composite absorption by 
multiple structures in a single line of underlines the diffi- 
culty of using the kinematics of DLAs to infer their nature 
(Prochaska & Wolfe 1997; Haehnelt, Steinmetz & Rauch 
1998). 

3.4 Does a maximum likelihood estimation of 
absorber give an accurate DLA size? 

Given the clumpy distribution of DLA-column density gas 
in our simulated galaxies, it is interesting to ask whether 
statistics from pairs of QSOs can give an accurate picture of 
galaxy 'size'. We therefore prepare a mock catalog of pairs 
of sightlines from which we calculate a maximum likelihood 
estimate of absorber size based on the coincidences and anti- 
coincidences of DLA absorption (e.g. McGill 1990; Dinshaw 
et al. 1997). Whilst relatively simple to implement, this tech- 
nique assumes a simple geometry (usually a spherical halo or 
cyclinder) and does not account for density fluctuations and 
other small scale structure. Nonetheless, this technique has 
previously been used for C IV and Mg II absorbers (Lopez, 
Hagen & Reimers 2000; Ellison et al. 2004), but the statistics 
are not yet available to attempt this for DLAs. For pairs of 
sightlines intersecting spherical clouds, the probability that 
a halo is intersected by the second line of sight, given that 
it is seen in the first, is given by 

</>(X) = - jarccos [X{z)\ - X(z) y/l - A(z) 2 } (2) 

for ^ X(z) ^ 1 and zero otherwise. Here, X(z) = 
S(z)/2R where S(z) is the line of sight separation and R 
is the absorber radius. In order to make maximum use of 



the information available we actually want to calculate the 
probability that both lines of sight are intersected, if either 
line of sight shows absorption. This probability is given by 



And the likelihood function as a function of radius is 
given by 

C(R) = l[4> [*(%)] n {1 - i> [*(*,-)]} - (4) 

i i 

where i and j denote the number of coincidences and 
anti-coincidences. The success of this technique relies on 
having a sufficient number of coincidences as well as anti- 
coincidences to constrain the likelihood function. From a 
mock catalog of 50 (half each from the XZ and XY orienta- 
tions) pairs of sightlines with the impact parameter of the 
first DLA absorber < 30 kpc from the centre of the galaxy 
and with all line of sight separations < 30 kpc, we have 
5 coincidences and 45 anti-coincidences for galaxy K15 at 
z = 3.0 (middle left panels of Figures [9l and [10|| . The results 
of the maximum likelihood analysis for this mock catalog 
are shown in Figure [TBI we determine a most likely radius of 
16.8^ 1 '| /ifp 1 kpc (95% confidence limits). A smaller sample, 
or wider pair separations would have led to fewer coinci- 
dences and hence poorer constraints on the absorber size. 
A visual comparison of this most likely size with the spatial 
distribution of gas in Figure [9] shows that the maximum 
likelihood method gives a quite reasonable estimate of the 
extent of DLA gas even when the distribution is clumpy; the 
covering fraction of gas with log N(HI) 20.3 within 16.8 
kpc is ~ 50%. Of course, a real catalog of pairs will likely 
probe galaxies of very different sizes and will therefore yield 
an average value. For a smoother morphology, such as the 
regular disk formed by the K33 galaxy at z = 2.3 (see bot- 
tom right panel of Figure [9} we find a most likely radius of 
12. 9^2 6 kP c w ith 90% of the gas within this radius above 
the DLA N(HI) threshold (XY plane). 

The compilation of a sample of 50 pairs with separations 
< 30 kpc (~ 4 arcsecs) in which at least one DLA is detected 
is a considerable challenge. However, the required number of 
pairs is fewer if the separations are also smaller, because co- 
incidences will be detected at a higher rate and therefore 
constrain the absorber size better for a smaller sample. For 
example, if we repeat the maximum likelihood analysis de- 
scribed above for the K15 galaxy at z = 3.0, but with a 
maximum pair separation of 20 kpc, from 4 coincidences 
and 16 anti-coincidences we deduce a most likely radius of 
13.2±2° 5 5 kpc. This means that if the simulations we have 
used in this work are a reasonably accurate representation of 
the size of the DLA cross-section in high redshift galaxies, 
close pairs with separations < 3 arcseconds or wide sepa- 
ration lenses will be best able to observationally constrain 
DLA size. 



4 CONCLUSIONS 

We have presented moderate resolution spectra of a pair of 
QSOs at z ~ 3 separated by 13.8 arcseconds on the sky. DLA 
or sub-DLA absorption is identified at three intervening red- 
shifts: z abs = 2.47, 2.66 and 2.94. For the two higher redshift 
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Figure 16. Maximum likelihood distribution (solid line) and cu- 
mulative distribution (dashed line) of the size of DLA absorbing 
gas from a mock catalog of 50 pairs of sightlines for galaxy K15 
at z = 3.0. The most likely radius from this simulated catalog is 
found to be 16.8 kpc with 95% confidence limits of 15.4 and 
21.6 h'g kpc. 



systems, absorption is seen in both lines of sight which have 
proper transverse separations of ~ 110 h^g kpc. We mea- 
sure chemical abundances for three of the five absorbers de- 
tected in our ESI spectra and determine high metallicities 
for two of them: 1/5 and 1/3 Zq. Both absorbers also have 
quite high dust depletions and for the z a b s = 2.94 DLA we 
also find that there is no a element enhancement, [S/Zn] 
= +0.05, and measure a low fraction of molecular gas, log 
/(Hh) < —5.5 in this particular line of sight. 

Although this is the first time that transverse DLA ab- 
sorption on this scale has been studied, if line of sight blend- 
ing occurs as often as the coincidences observed in SDSS 
1116+4118 AB, we demonstrate that this could have impor- 
tant ramifications. For example, the determination of chemi- 
cal abundances can be affected by several tenths of a dex and 
kinematics become challenging to interpret. Although more 
QSO pairs with separations of a few tens to a few hundreds 
of kpc are clearly required to confirm whether blending may 
really occur frequently in single line of sight DLAs, our work 
identifies a potential complication in the interpretation of 
DLA properties. 

This is the first time that DLAs have been studied on 
a transverse scale > 10 kpc and offers an opportunity 
to study their size and environment. By producing artifical 
pairs of sightlines through high resolution galaxy simula- 
tions, we conclude that the coincident absorption is unlikely 
to be associated with a single large galaxy, or a galaxy plus 
satellite system. Instead, from a clustering analysis, we find 
that the statistically more likely scenario is one in which 
the binary sightlines intersect multiple galaxies. However, 
the probability for coincident sub-DLA absorption is still 
< 10%, so more observations of QSO pairs are required to 
more robustly interpret our data. 

The prospects for extending the analysis of absorption 
in QSO pairs is promising, thanks to the large sky coverage 
of the SDSS. Hennawi et al. (2006b) estimate that the cur- 
rent 8000 deg 2 of SDSS imaging contains ~ 140 pairs with 
A6 < 25 arcsecs and z > 1.8. Many of these will require 



follow-up spectroscopy since the blue coverage of the SDSS 
spectra only reaches down to z a b s ~ 2.2. However, based 
on absorber number densities, Hennawi et al. (2006b) es- 
timate that the SDSS sample of pairs would contain ~ 70 
sub-DLAs or DLAs. Measuring the rate of coincident ab- 
sorption in a sample of this size may, at last, place the scale 
of DLA absorption and the structure of its environment on 
solid ground. 
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